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Abstract 

We apply effective field theory and renormalization group techniques to the problem 
of Cooper pair formation in neutron stars. Simple analytical expressions for the ^5*0 
condensate are derived which are free of nuclear potential model dependencies. The 
condensate is evaluated using phase shift data from neutron-neutron scattering. 
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Neutron superfluidity plays an important role in the physics of neutron stars, affecting 
the neutrino cooling rate and heat capacity as well as starquake phenomena . However, the 
computation of the neutron-neutron (NN) condensate is a difficult problem and has resulted 
in a wide range of results p|, ^ ^, . The main reason is the exponential sensitivity of the 
condensate to the effective interaction, which is typically extracted in a model-dependent 
way from NN scattering data, and also depends on medium effects. 

In this paper, we use the renormalization group (RG) to construct the effective field theory 
of neutrons near the Fermi surface. In this description the size of the condensate is related 
to the position of a Landau pole in a running coupling near the Fermi surface (FS) 0, |]. All 
interactions other than the Cooper pairing interaction (a four-neutron operator restricted 
to kinematic points describing scattering of neutrons with equal and opposite momenta) 
can be shown to be irrelevant at the FS. Due to this simplification the evolution of the 
Cooper pairing interaction, and hence the location of the Landau pole, can be determined 
analytically. Medium effects are reflected in the RG evolution of Fermi liquid parameters 
such as the effective mass and four-neutron operator, whose initial form G{p^) is obtained 
by matching to NN phase shifts via an exact expression for the in-medium NN amplitude. 
We focus on the ^Sq condensate, although our techniques can also be applied to study 
condensates of higher angular momentum. 

We begin by reviewing the effective field theory description of Fermi liquids [0, ^. In this 
description we make a guess as to the form of the effective theory close to the Fermi surface. 
The obvious guess based on the dynamics of non-relativistic systems is that the theory is one 
of weakly interacting fermions: these are the dressed "quasi-particles" of solid state physics 
language. We will henceforth refer to these effective degrees of freedom as "neutrons", 
with the understanding that they could in principle be related to the bare neutrons in a 
complicated way. Rather than treating other degrees of freedom such as pions, deltas, etc., 
as propagating degrees of freedom we will integrate them out leaving a potentially infinite 
sum over non-local, higher dimension fermion operators. The effective Lagrangian is simply 

C = i)i{i^ + /i7o - m)ipi + ■■■ , (1) 

where the ellipsis denote higher dimension interaction terms. Although we will eventually 
specialize to the non-relativistic limit, we begin with a relativistic formulation because it 
allows us to systematically track the corrections to that limit. 

The chemical potential in (0) naturally breaks the 0(3,1) invariance of space-time to 
0(3) and furthermore picks out momenta of order pp = y/Jj^^^-m^- It is therefore natural to 
study the theory as we approach the Fermi surface in a Wilsonian sense. We parameterize 
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four momenta in the following fashion 

p''=ip',p} = ik',k + l) (2) 

where k lies on the Fermi surface {\k\ = pp) and / is orthogonal to it. We study the Wilsonian 
effective theory of modes near the Fermi surface, with energy and momenta 

\ko\, |/1 < A , A . (3) 

While this type of RG scaling is somewhat unfamiliar, it actually corresponds to thinning out 
fermionic degrees of freedom according to their eigenvalues under the operator + /i7o — m. 
It is easy to see that eigenspinors of this operator with eigenvalues : |An| < A correspond 
to states satisfying (^. Consider an eigenspinor of the form 

= uiE,p} e'(P-P<^'^ , (4) 

where u{E,p} satisfies the usual momentum-space Dirac equation with E = ^p^ + = 
m + fi ± ^(^)) ^'^d pq = m ± 0{A). Then by direct substitution we see that tpp satisfies 

(z^ + /i7o - m)ilJk = 0(A) . (5) 

Thus, the RG flow towards the Fermi surface just corresponds to taking the cutoff on eigen- 
values of + /i7o — m to zero. 

Which operators are relevant in this limit? For our guess to make sense the kinetic term 
for the fermions must be invariant when we scale energies and momenta, ko — > sko (or, 
t —>■ t/s), and I — > si, with s < 1. We must be careful to satisfy all the global symmetries 
of the theory. In particular, there is a spurious symmetry of (|l|) in which we treat n as 
the temporal component of a fictitious U{1)b gauge boson. In other words, the combined 
transformations ip e*^* ip and fi ^ fi + 6 leave the lagrangian (|I|) invariant. From this, 
we deduce that time derivatives acting on ip and factors of fi may only enter the effective 
theory in the combination {Oq + /i)7o. This requires the kinetic term of our effective theory 
to be of the form 

'^e// = / '^^ ^ ((*^* + '")^o - p-i - rri^ip , (6) 

where m is the effective neutron mass. In the Wilsonian RG scaling, we eliminate all modes 
with energy and momenta |/i;o|, |^| > A, where A is our cutoff. As discussed above, on the 
remaining degrees of freedom the operator {idt + fi)'~fo—p-^ — m ~ 0{A) and therefore scales 
like s. We deduce that for (^ to remain invariant, -ip must scale as s"^/^. Now consider the 
four fermion operator 

G [ dt d^pi d^p2 d^ps d^pi ^F^ i)Ti) S\pi +P2+n+ Pi) (7) 
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where T contains any Lorentz or flavor structure. Naively, for / close to zero, the delta 
function does not scale, and is irrelevant since it scales as s. Higher dimension operators 
with extra powers of the fields are clearly irrelevant as well. The only operators that survive 
are four-fernaion operators satisfying the kinematic constraint pi ~ ~P2, in which case the 
delta function becomes 

S{h + l2-h- k) (8) 

and scales as s^^. The resulting four fermion operator is marginal, and quantum effects must 
be considered to determined its relevance to dynamics at the FS. 

We now study the RG flow of the marginal Cooper pairing operator. We are primarily 
interested in the ^5*0 component of this operator, which will be seen to dominate the others 
near the FS. We first consider the ^5*0 component by itself, and later consider corrections to 
this approximation that result from loops involving higher angular momentum components. 
It is easy to show that angular momentum conservation forbids such corrections so the ^5*0 
evolution is exact. 

Let us briefly review the kinematics of neutron-neutron scattering at the Fermi surface in 
the center of mass frame. Let the incoming momentum be pi, and the outgoing momentum 
f)2-, both of magnitude pp. The scattering amplitude is in general a function of the angle 
between these two vectors. The form factor of an operator describing this scattering process 
can be decomposed into components corresponding to projections onto different angular 
momentum eigenstates, or spherical harmonics y^(6',0). Actually, since the amplitudes are 
independent of only the m = components are required, and the projection just involves 
integration over 6. The ^5*0 component is obviously the component which is independent of 
6', and hence has a constant form factor. 
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The beta function for our four-neutron operator can be deduced from the one loop graph 
in figure If we consider only the ^Sq component of the coupling the powers of G can be 
factored out of the integral, yieldingF] 



We can rewrite this as 



ik 



(9) 



G' 



(/i7o + m)ifc (/i7o + m)j7 1 2 n / n 

^ ^ + Y^PF2^ila)ik{la)jl 



a=l 



(10) 



where the integral / is given by 



I = ^ 



Pf 



dQk f dkodl 



(27r)3 
A 



27T kl 



In 
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Here A/^j and kyy are the infrared and ultraviolet limits of integration. In the usual Wilso- 
nian sense, the effects of modes in the shell between these cutoffs is summarized in the 
evolution of the coupling G. In the non-relativistic limit, where ^ — 0, (0) becomes 



bG = G^ 



7o + 1 



ik 



7o + 1 



3I 



m pp 
27r2 



t 



(12) 



where t = ln(-^^). To incorporate non-relativistic corrections to some order in pp/m'^, one 
must include additional operators appearing in (^) in the RG equations. However, in our 
case of interest, p'jp/m'^ is at most of order a few percent, so we will drop all non-relativistic 
corrections. Since ^(70 ± 1) is simply the neutron/antineutron projection operator, the 
result only renormalizes the coupling between neutrons, and does not involve anti-neutrons. 
Henceforth we will focus on this interaction. The resulting RG equation is 



dG 



m Pp 
27r2 



G^ 



The solution to this equation is 



G{t) 



27r2G(0) 



which has a Landau pole at 



t 



27r2 - 1 mppG{0) 



mppG{0) 



(13) 



(14) 



(15) 



^Note that the ie prescription for propagators at finite chemical potential is slightly subtle Q. Essentially, 
the sign of the ie terms is such that the usual Wick rotation can be made regardless of the sign of {E{p) — /i). 



5 



By dimensional analysis, we expect a Cooper pairing gap of size 

2^2 

A ^ Xuve^^^^^ . (16) 

The exponent in ([T6|) is very similar to the usual BCS weak-coupling result 0, except that 
in our case the coupling G has a well-defined origin: it arises from the matching of our purely 
nucleonic effective theory to the full theory and can be extracted from NN scattering data. 
The prefactor will be determined by matching to known results for the gap at low density. 

Now we return to the issue of higher angular momentum components. As discussed 
previously we work in the basis provided by the spherical harmonics y^(6', 0). Actually we 
only require the m = harmonics, which are independent. Breaking G{6) into its spherical 
harmonic components 

G{e) = 0^"^ Y^{e,o) + G^^^Y^\e,o) + G^^^Y^{e,o) + ■■■ (i7) 

and repeating the previous analysis, the only change is in the replacement 

G2 Jd'k^ J (fk G{P,k)G{h,) (18) 

where (5ah is the angle between a and h. 

Inserting (|l^ and simplifying yields the new set of RG equations 



m 



..^p,Y.iG^'^fx'' (19) 

where x^'' is given in terms of integrals over Legendre polynomials. It is easy to show that 
x"' is diagonal, using the following result: 

t/2A;Fj(/?pfc,0) Fj'(/3fcg,0) (X rj(/3p„0) 5"' . (20) 

This implies that there is no mixing between different angular momentum components, and 
our treatment of the ^Sq RGE is exact. 

The final step in our analysis is to match the effective neutron Lagrangian to the full 
nuclear theory, which must include pions as well as neutron-neutron contact terms |T0 



However, since our FS effective theory includes only very low energy quasiparticles, the 
pions have been integrated out. We need only retain the non-local four- neutron operator, 
whose value can be fixed by comparison with ^5*0 NN scattering data. 

First we note that the ^5*0 scattering amplitude in the theory with only neutrons can be 
found by summing a bubble chain of Feynman graphs (figure ^ in which the vertices are 
given by the Lagrangian form factor G(p^, v). Here we work in the center of mass frame, so 
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X + XX+X)CX+ 



Figure 2: Bubble Sum leading to iJ^ 

G is only a function of the neutron momentum p. We perform our calculation in an effective 
theory living in a thin shell around the Fermi surface, regulated by a hard ultraviolet cutoff 
Auv and a hard infrared cutoff Am. The resulting amplitude can be directly computed and 
is given by 

1 + m z/) (i/ + i-p) 1 4:71 ' 

where 

^=—H-. — ) • (22) 

TC Auv 

Invariance of the physical scattering amplitude Ai under changes in u implies the following 
RGE for the coupling z/): 

1 1 miu — v') , . 

+ ^1 , (23) 



G(p2,z/') G{p^,u) Att 
and a corresponding Landau pole at 

Comparison with ( p!3[ ) reveals that the two RGEs are identical. The scale associated with 
the Landau pole is 

A^^ = kuve'''"'^'' . (25) 

Note that, unlike the scattering amplitude A^, the coupling GijP', z/) is not a physical quantity 
and cannot be determined without fixing a subtraction scheme (i.e. specifying v). This leads 
to an ambiguity in the overall normalization of the scale associated with the Landau pole. 

We can rewrite (plD in terms of a phase shift, defined in terms of the S-matrix hy S = e^^^: 

5 = \: XnU+i"^ m\ . (26) 



2i V 27r 

Finally, we can invert this relationship to obtain an expression for the coupling, 

G{p^,v) = '- — — . 27 

' ' u + p cot(5) ^ ^ 
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Substituting this into ( ^5]) yields an equation relating the superfluid gap (which must be 
equal to the scale of the Landau pole, up to a some factor of order one) to the phase shift, 

A = Ae-^'^"*^ . (28) 

The constant A is undetermined due to the ambiguity mentioned above, but scales like the 
Fermi energy ep = p|^/2m*, since it is proportional to the UV cutoff or FS shell thickness. 
The precise numerical value of the coefficient can be determined by studying the weak 
coupling limit of a low-density neutron gas [pp —>■ 0). Once determined, the coefficient 
remains fixed independent of the Fermi momentum pp. An explicit computation using the 
gap equation has been performed in this limit by Khodel et al. 0], with the result A = ^ej?. 
The final result is 

A = ^ Ae-t™*^ . (29) 
e2 2m, ^ ^ 

Our calculation can be viewed as a rigorous justification for the use of a BCS-like expression 

depending on Fermi liquid parameters which incorporate in-medium effects [0. In our case 



the coefficent of the exponential can be justified from first principles. Equation (^91) is exact 
when evaluated using the effective mass and in-medium phase shift data. If such data is not 
available, an estimate can be obtained by using the corresponding vacuum data, although 
as is shown in figure |^ the difference could be significant. 

In figure ^ we display the resulting g function of Fermi momentum pp, obtained by 

inserting ^So phase shifts obtained from the Nijmegen partial wave analysis of NN scattering 
data [jlll. The result of (|29|), using the bare neutron mass and zero-density (non- medium) 



phase shifts, is given by the solid curve. Inserting the effective mass and scattering ampli- 
tudes from [0 yields the short-dashed curve. The long-dashed and dot-dashed curves give 
the results of @] (lower curve) and (upper curve). 

It appears that the result (|29|) is accurate as well as simple. Our results agree reasonably 
well with state of the art calculations using more traditional methods. Our leading result 
(ignoring medium corrections) is intermediate between calculations with "induced potential" 
effects 1^, 0, 0, which have a maximum gap size of order 1 MeV, and those using a bare 
potential 0, |^, |^ , which produce a maximum gap of approximately 3 MeV. When medium 
effects are included the gap is decreased significantly. 
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Figure 3: Delta vs. Fermi momentum pp. All units are MeV. Solid curve represents Eq. 



(|29D. Long-dashed and dot-dashed curves are from and 0. Short-dashed curve is 



with inputs from 
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